Crossover dark soliton dynamics in ultracold one-dimensional Bose gases 
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Ultracold confined one- dimensional atomic gases are predicted to support dark soliton solutions 
arising from a nonlinear Schrodinger equation of suitable nonlinearity. In weakly-interacting (high 
density) gases, the nonlinearity is cubic, whereas an approximate model for describing the behaviour 
of strongly-interacting (low density) gases is one characterized by a quintic nonlinearity. We use 
an approximate analytical expression for the form of the nonlinearity in the intermediate regimes 
to show that, near the crossover between the two differ ent regimes, the soliton is predicted and 
numerically confirmed to oscillate at a frequency of y2/3fi, where D. is the harmonic trap frequency. 



I. INTRODUCTION 

Dark solitons (DSs), the most fundamental nonlin- 
ear excitations of the one-dimensional defocusing non- 
linear Schrodinger (NLS) equation, have been studied in 
a broad range of physical systems. Apart from the the- 
oretical work, experimental studies on DSs include their 
observation either as temporal pulses in optical fibers [ll , 
or as spatial structures in bulk media and waveguides [2j 
(see also Q for a review), the excitation of a nonpropa- 
gating kink in a parametrically-driven shallow liquid [j] , 
standing DSs in a discrete mechanical system [5|, high- 
frequency DSs in thin magnetic films [6|, and so on. 

Recently, DSs have attracted much attention in the 
physics of atomic Bose- Einstein condensates (BECs) [7| , 
where dark matter-wave solitons have also been observed 
experimentally 8]. Dark solitons in BECs arc known 
to be more robust in one-dimensional (ID) geometries 
and at very low temperatures (a regime which is cur- 
rently experimentally accessible @). For that reason, 
the majority of theoretical studies on DSs in BECs have 
for simplicity been performed in the framework of the 
one-dimensional (ID) cubic NLS equation, which in this 
context is referred to as the Gross-Pitaevskii equation; 
the latter, is the commonly adopted mean-field theoretic 
model describing ultracold weakly-interacting Bose gases 
in the absence of thermal or quantum fluctuations. Many 
of the above mentioned theoretical studies have been de- 
voted to the analysis of the dynamical properties of mov- 
ing DS, such as their oscillations [l(| El, El, El, EH and 
sound emission [ll], El, EH in the presence of the ex- 
ternal trapping potential. On the other hand, and in 
the same context of the ID Bose systems, DSs have also 
been studied in the framework of a quintic NLS equa- 
tion 

BE! El llTjl. a long- wavelength model which has 
been proposed [1 51 ] for the opposite limit of strong in- 
teratomic coupling [IH [l9j ; in this case, the collisional 
properties of the bosonic atoms are significantly modified, 
with the interacting bosonic gas behaving like a system 
of free fermions [2fJ. Such, so-called, Tonks-Girardeau 
(TG) gases have recently been observed experimentally 
as well [HI. 



Regarding the dynamical features of the moving DSs 
in trapped ID Bose gases, the above works revealed that, 
in the absence of other dissipative mechanisms, the soli- 
ton oscillates in the trap with a frequency which differs 
between the weakly and strongly interacting regimes. In 
particular, in the presence of a harmonic confining poten- 
tial of frequency fi, the study of the cubic (quintic) NLS 
predicts such an oscillation frequency to be afl. where 
a = 1/^/2 (a = 1) for the weakly 0, El El, El El 
(strongly [3 EH) 

interacting Bose gas. The latter result 
for the strongly interacting case is identical to the corre- 
sponding one obtained by a full many-body calculation 

The transition between the weakly and strongly inter- 
acting regimes is usually characterized by a single param- 
eter, denoted by 7, quantifying the ratio of the average 
interaction energy to the kinetic energy calculated with 
mean field theory [19( . This parameter varies smoothly as 
the interatomic coupling is increased from values 7 <C 1 
(weakly interacting regime), to 7 ^> 1 (strongly inter- 
acting regime); thus, an approximate "crossover regime" 
can be identified around 7 ~ 0(1), as also attained ex- 
perimentally [23|] ■ The value of 7 can be controlled ex- 
perimentally by various independent parameters, such 
as scattering length, transverse confinement, density, or 
even modification of the effective mass of the system. 

Motivated by the investigation of the DS dynamics 
in the two limiting regimes, in the present work we ob- 
tain analytical results, which are confirmed by numerical 
simulations, for the oscillation frequency of a DS in the 
crossover regime for a purely ID system. In particular, 
we show that for 7 ~ 0(1) the DS oscillates with a fre- 
quency which is higher than the one (-\/l/257) 
pertaining to the weakly interacting Bose gas. Such an 
increase in the soliton oscillation frequency could serve as 
an additional diagnostic test for the deviation from pure 
bosonic mean field behaviour. 

The paper is organized as follows: In section II we 
present the generalized NLS model proper and find its 
ground state, as well as its linear (sound waves) and 
nonlinear (dark solitons) excitations. Section III is de- 
voted to the analytical derivation of the DS oscillation 
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frequency and the discussion of relevant numerical simu- 
lations. Finally, in section IV we summarize our findings. 



II. THE MODEL AND ITS ANALYTICAL 
CONSIDERATION 

A. The generalized NLS equation 

Importantly for our present analysis, the parameter 7 
separating the regimes of weak and strong interatomic 
coupling can be re-expressed, for a given system configu- 
ration, in terms of the inverse ratio of the system density 
to some critical density. This enables us to perform a uni- 
fied analysis of all regimes by means of a NLS equation 
with a generalized nonlinearity, for the parameter ip con- 
nected to the density n = \ip\ 2 of an ultracold confined 
ID Bose gas. This equation takes the form 



..dip 



2m dx 2 



+ V ext (x) +$(n) 



(i) 



where m is the atomic mass and V ey ±(x) — (l/2)muj 2 x 2 is 
the external trapping potential (u) x being the axial con- 
fining frequency). 

The exact form of the nonlinearity <I>(n) valid in both 
limiting regimes and the crossover region is well-known in 
the homogeneous hydrodynamic limit. While the func- 
tional dependence of $(n) on 7 (and its analytical asymp- 
totics) are known [3] , its precise values in the "crossover 
region" can only be evaluated numerically. Such inter- 
mediate values have been tabulated by Dunjko et al. in 
[HI], and subsequently discussed by various authors in 
the local density approximation, see e.g. [H, |2f|. Since 
we are interested in this "crossover region" , we should 
thus use an approximate expression for the nonlinearity, 
which captures both limiting regimes exactly and pro- 
vides a good approximation for intermediate values. At 
the same time, however, we are constrained in our present 
work by the need for a relatively simple expression which 
will enable rather involved analytical work to be carried 
out. For our purpose, it is thus sufficient to use a some- 
what simplified generalized nonlinearity <3?(n) of the form 



$(n) 



2m 1 + n/h c 



(2) 



which nonetheless corresponds to a fairly good analytical 
approximation [26|, [27j to the exact nonlinearity [HI . In 
this notation, the critical density approximately mark- 
ing the crossover region is given by n c — 4/(7r 2 an)), 
where am is the effective ID "scattering length". This 
is defined in terms of the usual three-dimensional s- 
wave scattering length 0,30 via cl\d = l\/a>&Di where 
l±_ = ^Jh/moj± is the harmonic oscillator length in the 
transverse direction. 

Equation |T]) is well-documented in the high-density 
limit n^> n c (corresponding to 7 <C 1). In this case, the 



parameter ip describes the "wavefunction" of a weakly- 
interacting ID Bose gas under harmonic confinement, for 
which Eq. (TTJ) reduces to the ID Gross-Pitaevskii equa- 
tion, 

dip h 2 d 2 TA , Airti 2 ( a 3D 

2m dx 2 



ih- 



dt 



(3) 

In the opposite limit of strong interatomic coupling, 
which, rather counter-intuitively, corresponds to the low 
density limit n <C n c , the parameter ip satisfies the fol- 
lowing quintic NLS equation, 



ih 



dip 
"dt 



h 2 d 2 

2m dx 2 



+ V cxt (x) + 



.2fc2 



2m 



^ (4) 



which has been derived by various different methods in 
[H HI]. While the validity of Eq. (jH) to discuss co- 
herence properties of strongly-interacting ID Bose gases 
has been questioned (29j . the corresponding hydrody- 
namic equations for the density n and the phase (p (or 
the atomic velocity v = d x (p) arising from this equa- 
tion under the Madelung transformation ip = v /nexp(i0) 
are well-documented in the context of the local density 
approximation [25[. An equation of the form Q, 
which explictly includes the so-called "quantum pressure 
term" (H 2 /2m)d 2 ip /dx 2 , should however only be valid for 
density variations which occur on a lengthscale which is 
larger than the Fermi healing length £ = l/(7m p ), where 
rip = I ip (0)| 2 is the peak density of the gas at the trap 
center. To avoid such potential complications, the physi- 
cal analysis presented in this paper is only concerned with 
the limit of shallow DSs, for which there is only a very 
slow density variation within the Fermi healing length. 

Measuring the variables x and t, and the density \tp\ 2 , 
in units of the Fermi healing length £, the time tq = 
h/mt; 2 , and the density v2n p respectively, Eq. fl} can 
be rewritten in the following dimensionless form, 



iu t = ~^u xx + V(x)u + g(n)u, 



(•5) 



where the subscripts denote partial differentiation. In 
view of the above scalings, the normalized confining po- 
tential becomes V{x) = (l/2)(£ i /l x ) 4: x 2 , where l x = 
yhjmujx is the harmonic oscillator length in the ax- 
ial direction. As the parameter £/l x is apparently 
small, it is convenient to define the small parameter 
e = Sl~ 2 / 3 (£/Z) 4 / 3 where fl is a parameter of order O(l), 
that will be used in the perturbation analysis to fol- 
low. This way, the external potential is actually a func- 
tion of the slow variable X = e 3 ^ 2 x, and has the form 
V(X) — (1/2)Q 2 X 2 , where expresses the trap fre- 
quency. Finally, the nonlinearity function g(n) (with 
n = \u\ 2 being the normalized density) becomes, 

t \ n2 u 2V2 1 

9\ n ) = 1 I TZ » where n c = — — , (6) 



l + n/n c ' tt 2 aiorip 

The weakly and strongly interacting limits discussed 
above respectively arise when the dimensionless param- 
eter 7 = 2/ {nam) [3 obeys 7 < 1 or 7 > 1. In our 
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present notation, n/n c = (tt 2 /2)j~ 1 , so that our approx- 
imate crossover region actually corresponds to a value 
7 = 2/ir 2 w 0.2. For such a relatively small value of 7, it 
is still reasonable to use Eqs. ©-© [or Eqs. ©-©] to 
describe this physical system. 

Studies of DSs based on generalized ID NLS equations 
first appeared in the literature for homogeneous systems 
to deal with saturable nonlinearities appearing in the 
context of nonlinear optics in this case, the nonlin- 
earity becomes cubic at low densities, rather than at high 
densities which occurs for ultracold pure ID Bose gases. 
More recently, DSs in generalized ID NLS equations were 
also considered in the BEC context, but only as effec- 
tive theories for weakly-interacting elongated 3D conden- 
sates. These equations contained either a cubic-quintic 
nonlinearity with constant coefficients (30j . or a gener- 
alized non-polynomial nonlinearity depending explicitly 
on the trap aspect ratio [3l|, H3] ■ The aim of this work is 
rather distinct, namely to calculate the soliton oscillation 
frequency in the crossover between weakly- and strongly 
interacting ID Bose gases. The validity of our work is 
thus restricted to the pure ID regime, with the precise 
experimental conditions needed for such a crossover dis- 
cussed in detail in [33| . 

Below we apply the reductive perturbation method 
(see, e.g., [34L l35|). valid in the limit of shallow soli- 
tons, to analytically obtain the soliton oscillation fre- 
quency close to the critical density n c , which approxi- 
mately marks the crossover between the regimes of weak 
and strong interatomic coupling. 



B. Ground state, linear and nonlinear excitations 

Starting off from the generalized NLS Eq. ([5]), we use 
the Madelung transformation, u = -Jn exp(i(ft), to obtain 
the following set of hydrodynamic equations, 



n t + {n(ft x ) x = 0, 
<k + 9(n) + -4>l 



(7) 

^n- 1 / 2 (n 1 / 2 ) xx + V(X)^0.(8) 



The above equations are similar to the ones that have 
been employed to discuss the crossover from TG to BEC 
regime 18, 25]. The ground state of the system can 
be obtained upon assuming that the atomic velocity 
v = 4>x — (i.e., no flow in the system) and (ftt = — /io 
(dimensionless chemical potential). Then, as Eq. ([7]) im- 
plies that n = tiq is time-independent in the ground state, 
we assume that no — uq(X). Thus, to leading order in e 
[i.e., to 0(1)] Eq. © yields, 



g(n ) = fio- V(X), 



(9) 



in the region where /io > V(X) and hq = out- 
side. Equation ([9]) determines the density profile in 
the so-called Thomas-Fermi (TF) approximation; note 
that for the typical case of the harmonic trap, e.g., 
V(X) = (l/2)fl 2 X 2 , Eq. © recovers the well-known 



result that the density profile is parabolic, with iio(X) = 
n~ x [no — V(X)], in the regime n> n c and elliptic, with 
no(X) = \J /io — V(X), in the regime n <C n c . Also, it 
is noticed that from Eq. ([9]), and for the harmonic trap 
under consideration, the axial size of the gas is 2L, where 
L = ^/2/i /f2 is the TF radius. 

We now consider the propagation of small-amplitude 
linear excitations (e.g., sound waves) of the ground state, 
by seeking solutions of Eqs. ©-© of the form n = 
no(X) + en(x,t) and (ft = —jJ,ot + €(ft{x,t), where the 
functions h and (ft describe the linear excitations. In- 
serting this ansatz into Eqs. 0-ldJ), to order 0(1) we 
recover the TF approximation, while to order O(e) we 
obtain a system of linear equations for the linear excita- 
tions. Assuming plane wave solutions of this system, i.e., 
(n, (ft) ~ exp[i(kx — u>t)], we readily obtain the dispersion 
relation lu 2 = gonok 2 + fc 4 /4, where go = (dg / dn)\ n=no . 
This dispersion relation has the form of a Bogoliubov- 
type excitation spectrum, but with the excitation fre- 
quency u> being a function of the slow variable X. The 
speed of sound is local, due to the presence of the external 
potential, and is given by 



C = \fgan Q , 



(10) 



Note that Eq. (|10p shows that the speed of sound is 
given by C — y/n c riQ for n 3> n c , and C = V2n>o in the 
opposite regime n <§; n c . 

Next we analyze the evolution of the nonlinear excita- 
tions on top of the ground state, employing the reductive 
perturbation method [34], [35[ (see also [Tl| and [l6| for 
relevant studies in Bose gases). As the system of Eqs. 
©-([H]) is inhomogeneous, we introduce a new slow time- 
variable T = e 1 / 2 (t ~ f£ C^ix^dx'), and the following 
asymptotic expansions for the density n and phase (ft, 



= n {X) + en 1 (X,T) + e 2 n 2 {X,T) + 
= -/! < + e 1/2 <MA,T) + 



e 3/2 MX,T) 



■ -(H) 

Substituting the expansions (fTT|) into Eqs. ©-©, we 
obtain the following results: First, to order O(l), Eq. 
© leads to the TF approximation [see Eq. ©]. Then, 
to the first-order of approximation in e [i.e., to orders 
0(e) and 0(e 3/2 )], Eqs. © and © yield the equation, 



h(X,T) = -g (X) [ n x {X,T')dT', 
Jo 



(12) 



connecting the unknown functions n\ and (ft\. Finally, to 
the second order of approximation [to order 0(e 2 ) and 
0(e 5 / 2 )], Eqs. © and © lead to the following nonlinear 
evolution equation for rii, 



nix 



(3g + n g ) n 
2C 3 

ln(|C|3o) 1/2 



mniT 



1 

8C5 



(13) 



d 

~ ~dX 

where ijo = {d 2 g / dn 2 )\ n=nQ . Equation (fT3")) has the form 
of a Korteweg-deVries (KdV) equation with variable co- 
efficients, which has been used in the past to describe 
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shallow water-waves over variable depth, or ion-acoustic 
solitons in inhomogeneous plasmas [351 ]. Moreover, such 
KdV equations have recently been used to analyze the 
dynamics of DSs in Bose gases both in the weakly- 
interacting (llj and the strongly-interacting [l6| regimes. 

As the inhomogeneity-induced dynamics of the KdV 
solitons has been studied analytically in the past [36[ , we 
may employ these results to analyze the coherent evolu- 
tion of DSs in the Bose gas under consideration. Thus, 
introducing the transformations \ — f(8C) dX and 
ni = (3/2) (3(?o + n Q g y 1 C- 2 v( X ,T), we first put Eq. 
(fl~3l) into the form, 



v x — Qvvt + vttt = A(x)f, 



(14) 



where A(x) = (d/dx) In 



3/4 A l/4 
l 



In the 



»■>' 9o (3go + n g ) 
case A = 0, i.e., for a homogeneous gas with tiq(X) = 
n p = const., Eq. (fT4|) is the completely integrable KdV 
equation, which possesses a single-soliton solution of the 
following form [371 ] . 



v = -2K 2 sech 2 Z, Z = k[T — C(x)] 



(15) 



where £(%) = 4k 2 % + Co is the soliton center (with 
d(/dx = 4:K 2 being the soliton velocity in the T~x ref- 
erence frame), while k and Co are arbitrary constants 
presenting the soliton's amplitude (as well as inverse tem- 
poral width) and initial position respectively. Equation 
(|15jl describes a density notch on the backround density 
n p , with a phase jump across it [see Eq. (fT2")l . which 
implies that 4>i ~ tanh Z] and, thus, it represents an ap- 
proximate DS solution of Eq. ([5|). 

On the other hand, in the general case of the inho- 
mogeneous gas [i.e., in the presence of y(A)], soliton 
dynamics can still be studied analytically, provided that 
the right-hand side of Eq. (114[) can be treated as a small 
perturbation. As A(x) is apparently proportional to the 
density gradient, such a perturbative study is relevant 
in regions of small density gradients (e.g., near the trap 
center for a harmonic trapping potential), which is con- 
sistent with the use of the local density approximation. 
In this case, employing the adiabatic perturbation theory 
for solitons (38|, we may seek for the soliton solution of 
Eq. |H]) in the form of Eq. (jT5j), but with the soliton 
parameters k and £ being now unknown functions of x- 
The respective evolution equations for the soliton's am- 
plitude and center can be solved analytically [361 ] and the 
results, expressed in terms of the slow variable X, read: 



k(X) = k(Q) 



x 



2/3 



(16) 



1 



2k(0) 



1 - 



a(0) 



-2/3' 



(17) 



where k(0) and cr(0) are the values of the respective func- 
tions at X = 0. 



Note that the above procedure is general, and does 
not rely on the ratio of the parameter n/n c , although it 
has been shown to yield the correct results in both lim- 
its n 3> n c (ll| and n -C n c fl6j - In this work, we use 
the above general results for the evolution of the soli- 
ton parameters, to derive the equation of motion of the 
DS and find its oscillation frequency in the "nonlinearity 
crossover regime" n ~ n c . 



III. SOLITON OSCILLATION FREQUENCY 

Confining ourselves to the "crossover regime" n « n c , 
we may use a Taylor expansion of the function g(n) 
around n = n c , namely g(n) ~ g(n c )+g(n c )(n—n c ) + - ■ ■ , 
to derive the approximate expression go w n c (3no— n c ) / 4; 
the latter, along with the TF approximation [see Eq. ©], 
leads to the following density profile, 



n 



■■id 



4 \ L 



(18) 



where 5 = n c / \fjIo (it is reminded that L = y / 2/io/ri 
defines the axial size of the gas). Based on Eq. (fT8|) . it 
is now possible to derive the equation of motion of the 
DS as follows. First, we find the soliton phase, which, to 
order 0(e 3 / 2 ), reads, 



Z = e^ 2 K(X) 



t - 



dX' 
C(X') 



'^ 2 n(4 + s 2 ) 5/2 \^ 



(19) 

Then, looking along the characteristic lines of soliton mo- 
tion, it is possible to show that the position of the soliton 
satisfies the following equation of motion, 



dX 
~dt 



C 



1 



= (9\/3/2) K 2 (0)(4 At + S 2 )- 5 / 2 C' 



(20) 



For sufficiently small e the second term in the denomina- 
tor can be neglected; in this case, Eq. (|20|) shows that 
the velocity of sufficiently shallow DSs is approximately 
the same as the speed of sound given by Eq. lfTt)|) , i.e., 
dX/dt w \/gonQ. Thus, Eq. (|2H)) can be approximated 
by the separable differential equation, 



n 



-.dX = dt, 



n n c 



(21) 



which can readily be integrated. In particular, taking 
into account Eq. (fT5]) , we find that Eq. ([2"Tj) leads to the 
following result: 



4V13 
where h(S) 



h(S) \nw(X) + arcsin 
6(1 



n 



t, (22) 



¥ 2 ) 



1/2 



w{X) = 



x/l 2 - X 2 + ^h(S)X 

Vl 2 - x 2 



qzh(5)x 



(23) 
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and L = L^J 1 + |<5 2 . It can be seen that in the regions 

of small density gradient where X/L is sufficiently small 
[which is consistent with the assumption that the per- 
turbation ~ A(x) in the KdV Eq. (|14jl is weak], and for 
sufficiently small values of the parameter 5, the first term 
on the left-hand side of Eq. (|2"2")) can safely be neglected 
[39l ] . Then, it is readily seen that Eq. ([22]) is reduced to 
the following equation, 



X = L sin 




(24) 



Equation (|24|) demonstrates that a shallow DS displays 
an oscillatory motion in the harmonic trap V(X) = 
(l/2)n 2 X 2 in the "nonlinearity crossover regime" n f=a n c 
of Eq. (2), with an oscillation frequency Sl osc given by 



a 



(25) 



This finding has also been verified through appropri- 
ately crafted numerical experiments. As a typical ex- 
ample, we show, in particular, the evolution of a shal- 
low gray soliton (originally located at the origin) with 
initial speed C = 0.9 on the background of a potential 
V(x) = tl 2 x 2 /2, with n = 0.0707; n c was chosen to be 1 
in this case. Figure [T] shows the evolution of the space- 
time contour plot of the reduced density (the ground 
state density minus the actual density) for the NLS equa- 
tion, alongside our theoretical prediction of Eq. (|2"4")l . 
It is clear that during the first period, where the soli- 
tary wave does not significantly interact with the emit- 
ted "radiation" , the agreement between the theoretical 
prediction and the numerical results is very good. Sub- 
sequently, as expected given the above interaction, this 
agreement deteriorates. The approximate initial condi- 
tion is produced by imposing a NLS gray soliton on top 
of the ground state of the system for the nonlinearity 
of interest. Finally, we note in passing that our ana- 
lytical considerations are, strictly speaking, valid away 
from the turning points (where the wave speed vanishes). 
However, our numerical results, as well as alternative ap- 
proaches such as the ones of [l3[ and , illustrate that 
the range of validity of our results is, in fact, wider than 
what may be expected based on the mathematical limi- 
tations of the method. 

Note that the established values of the oscillation fre- 
quencies in the two limiting regimes of interatomic cou- 
pling can also be found in the framework of the presented 
analysis, upon utilizing Eq. (|2ip in the relevant limits 
n/n c 3> 1 and n/n c <C 1 and using the respective den- 
sity profiles. In the limit n ^> n c coresponding, e.g. to a 
weakly-interacting Bose-Einstein conde nsat e, the oscilla- 
tion frequency is O osc = £7/^/2 (or Vl a/2/4), whereas, in 
the opposite regime n < n c it has recently been found 
[3 O H1J that the oscillation frequency is £l osc = fl 
(or J7-\/2/2). The oscillation frequency in Eq. (|23|) is 
thus predicted to lie between the above mentioned lim- 
iting cases. This suggestes a continuous change in the 
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FIG. 1: Spatio-temporal evolution of the reduced condensate 
density (the ground state density minus the actual density) for 
the NLS equation with g(n) = n 2 /(l + n) and V(x) = S} 2 x 2 /2, 
with = 0.0707. The dashed line shows the theoretical pre- 
diction which agrees well with the full numerical result over 
(at least) the first oscillation period. 



predicted soliton oscillation frequency from the regime 
n > n c to n < n c . 

Assuming that the temperature is low enough for the 
soliton to perform many oscillations before decaying due 
to additional dissipative mechanisms excluded from the 
NLS equation, an observation of change in the oscillation 
frequency can be treated as a possible diagnostic tool of 
the system being in a particular interaction regime. A rel- 
evant experiment might thus be to create a Bose gas of a 
certain unknown interaction strength, and phase imprint 
a DS in such a system. Measurement of the oscillation 
frequency of this soliton will then provide important in- 
formation on the system parameter regime, with any de- 
viation from the oscillation frequency of denoting 
a regime of sufficiently strong correlations which, in turn, 
indicate that the weak-interaction model is no longer 
an adequate description of the system. Alternatively, 
one could create a dark soliton in a weakly-interacting 
ID BEC, and gradually increase the effective interaction 
strength, as done in some experiments. In this case, one 
should be able to observe a gradual monotonic increase 
of the oscillation frequency, even though the longitudinal 
confinement remains unaffected. 



IV. CONCLUSIONS 

In summary, we have discussed the dynamics of both 
linear and nonlinear excitations within a generalized non- 
linear Schrodinger equation motivated by considerations 
of the behaviour of ultracold atomic ID Bose gases. The 
considered model differs from relevant ones appearing in 
the context of optics |3j in that the linear behaviour in 
the density dominates at high densities, with quadratic 
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behaviour in the density dominating at lower densities. 
Within this generalised model, discussed here for the first 
time in relation to dark solitons in fD ultracold Bose 
gases, we have studied the dynamics of dark solitons in 
the "crossover regime" in the effective atomic interaction 
strength, which is approximately marked by a critical 
density. Our main conclusion, stemming from analyt- 
ical considerations and confirmed by numerics, is that 
the soliton oscillation frequency in this regime lies be- 
tween the known values arising in the two limiting cases 
of weakly and strongly interacting gases, indicating a 
continuous change between these two regimes. Finally, 
we note that although motivated by a particular physi- 
cal system, our model and analysis are quite general and 



are not restricted to the details of this particular system. 
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